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2H ' Masses of fermions in the SO(IO) 16-plet are constructed using only 

\^ . the 10, 120 and 126 scalar multiplets. The mass matrices are restricted 

CN I to be hermitian and the theory is constructed to have certain assumed 

quark masses, charged lepton masses and CKM matrix in accord with 
^^^ . data. The remaining free parameters are found by fitting to light 

t^ I neutrino masses and MSN matrices result as predictions. 

>'■ 

X 1 Introduction 

The simplest SO(IO) treatment of fermion masses, with the fermions being 
in SO(IO) 16-plets, uses the composition of 256 pairings of fermions into 
the 10, 120 and 126 SO(IO) represent at ions. There are 3 different 16-plets for 
the 3 generations of fermions and the fermion pairs are coupled to scalar 
(generalized Higgs) bosons in the same SO (10) representations with 3x3 
Yukawa coupling matrices; and when the scalar bosons develop vevs the 
3x3 fermion mass matrices are generated. The philosophy in this paper is to 
assume fermion masses only arise through coupling into these representations. 
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In the coupling to the {10} and {126} bosons SO(IO) gives equal coupling 
to the mass terms ipA'4'B and iPb4'a where A, B label different 16-plets of 
fermions. Here they are taken to be generation labels and thus the Yukawa 
coupling matrices can be taken to be symmetric without loss of generality; 
contrairiwise for the {120} the coupling of the mass terms ^aV's and V'bV'a 
are of equal magnitude but opposite signs so that the {120} Yukawas can be 
taken anti-symmetric [I]. For this, chiral, theory there is no constraint, in 
principle, that the Yukawa or mass matrix elements be real. 

In the very many SO(IO) fermion mass investigations, ranging from sim- 
ple hypotheses to complicated varieties of SUSY GUTS, there has been un- 
til recently a preponderance of hypotheses with no {120} Higgs particles 
[21 El in El E] • This rather arbitrary neglect seems to have been motivated, 
not unnaturally, by the search for simplicity, solvability and predictive power. 
Latterly, partly influenced by neutrino data, there have been more papers 
taking account of the {120} [HEIE] - some of these still giving a predomi- 
nant role to {10} and {126} ah initio. This present paper is an example of 
the {120} filling a necessary role, that of supplying the CP non- invariant, 
imaginary part, of the CKM matrix. 

There have been continual developments of SO(IO) type theories - among 
other reasons there have been the demands of data matching. Begining with 
theories having real Yukawa coefficients and real vacuum expectation values 
these have ranged to much more sophisticated models such as many having 
vevs with arbitrary phases to be determined from data extrapolated to a 
GUT scale (as for example in refs. [21 |8]). In this paper a simplification 
from general possibilities is chosen - that is the contribution to the fermion 
mass matrices from the {10} and {126} to be real and symmetric while that 
from the {120} be pure imaginary and anti-symmetric, resulting in Hermitian 
mass matrices. That condition on the {10} and {126} is familiar, especially 
in early papers, and can be presented as 'naturally' implemented using real 
Yukawas; in the formalism given in the next section the above condition 
on the {120} contribution can be presented as equally 'natural'; that is the 
mathematical formalism, as outlined in section 3, suggests the hermiticity as 
the simplest choice. 

There is also an important physical motivation for the choice of hermitic- 
ity, this being the association in the formalism with parity invariance at high 
energy. If the 75 terms are absent in the mass terms these are straightfor- 
wardly parity invariant, giving a limiting case of the chiral theory for the 
quark and charged lepton equations. And the absence of the 75 terms gives 



the hermiticity and then also the formahsm yields the mass contributions of 
the {10} and {126} as CP conserving and that of the {120} as CP violating. 
It can thus be said that there exists at least one basis (implicit in the formal- 
ism) in which the above mass matrices are hermitian. The one exception to 
high energy parity invariance and hermiticity in the mass terms of the model 
results indirectly, that is in the light neutrino masses from the contribution 
of the original seesaw (often known as 'Type V) in which the product of three 
hermitian mass matrices yields a non-hermitian matrix. The details for the 
model are in section IV below and the Appendix. (As has been emphasized 
recently [12] parity restoration in proceeding to high energies was one of the 
main motivations for the left-right symmetry principle which features in the 
next section.) 

The hermiticity, as derived below, is more precisely that of 3 x 3 flavour 
matrices associated with the symmetry representations, linear superpositions 
of which form the hermitian mass matrices of the quarks and charged leptons 
- and one other contribution to the light neutrino mass matrix 

At least two papers [TUl [H] incorporating the 120 along with hermitian 
mass matrices have previously been published. Both of these are set in MSSM 
theory and one in particular [10] has much detailed discussion of the various 
symmetry breakings in the MSSM context. The present paper is based on 
SO(IO) but has no committment to supersymmetry or details of a Higgs 
mechanism or other higher theory. 

The experimental data input to this model are the quark and charged 
lepton masses and the CKM matrix; so the parameters of the model have 
to be chosen to accomodate these numbers some of which carry considerable 
uncertainty; there can be further uncetainties from extrapolation to higher 
energies. The vital question then is can one adjust the very few remaining 
completely unknown parameters so that the model be compatible with the 
likely masses and MSN matrix of the three light neutrinos. In this paper 
the SO(IO) vevs (or substitute mechanism) giving rise to the masses are 
parameters to be chosen to match data and no particular Higgs potential (or 
other mechanism) is postulated.. 

In section II the SO (10) Clifford algebra formalism for the fermion masses 
is given. Section III outlines the calculation with the Hermitian matrix hy- 
pothesis and gives the resulting mass formulae in terms of Yukawa coefficients 
and scalar vacuum expectation values; while in section IV the assumptions 
of the theory allow the expression of the mass matrices in terms of Hermitian 
matrices and real ratios of vevs and also allow the incorporation of the quark 



and charged lepton mass data and CKM complex matrix data. Section V 
deals with the neutrino masses with emphasis on simple examples. One illus- 
trates tri-bimaximal mixing arising from neutrino masses suggested by the 
experimental data. 

2 SO(IO) and its 16-plets 

The Clifford algebra formalism of SO(IO) |1] is based on ten gamma matrices 
(Fi, ...., Fio) giving the 45 generators S^,^ = [F^, Vy\/2i. The gammas can also 
be expressed through the creation and annihilation operators XjiAjiiJ = 
1, ...5) where 

F2j_i = -i{xj - x]), ^2j = iXj + x]) 

The fermion 16-plet,?/'+ of positive lOd chirality is expressed as 
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|^+) = |0) ^0 + ^Xlxi |0) ^jk + -/""'"^-^xlxWrnXi |0) ^, 



where i/Jo, '^jk,'ipj are 2-component left-handed spinors of the particular gen- 
eration. The assignment to the leptons and the quark SU (3) colour triplets 
is, in an obvious notation with a, b being indices 1, 2, 3 and d, u being colour 
triplets: ^a ~ {dRY;ipa5 ~ d^, i^ab ~ {urY; ^a4 ~ ul; ^4 ~ ei;;V'45 ~ (e^)''; 
ip^ ^ m ;V'o ~ —{i^rY "where Ur denotes right-handed (heavy) neutrinos. 

The conjugate of equation ([T]) transforming appropriately under SO(IO) 
is 

(^:| Br ^ (0| i,^ + i,l (0| {\x]xh + V^J (0| Xj (2) 



where By = zFiFgFsFrFg and [^ = (0| XbXiX?.X2Xi- 

Let the suffices A and B be generation indices for the 3 generations. Then, 
letting C be the charge conjugation matrix 



{A\X\B)^{r+A BrC-'X\ij+B) (3) 

with X = F^ or X = F^F^^Fp or X = F^F^FpF^-F^ form the SO(IO) represen- 
tations 10 or 120 or 126 respectively. Under the action of a symmetry group 
generator T,^^ 

{A\X\B)^{A\[X,i:,,]\B). (4) 



Multiplying respectively by scalar fields 0^, 4>fj.up,(piiupaT of the same repre- 
sentations gives SO(IO) invariants, for example 



{A\ T, \B)<P^^ {r+A BrC-'T^ 1^+^) 0^ 

which for colourless neutral vevs contribute to the elements of the 3 by 3 
generation mass matrices of the quarks and leptons [1] . 

It is convenient to classify these vevs as the neutral colourless members 
of multiplets of the Pati-Salam subgroup of SO(IO). Selecting the generators 
formed by (ri,....,re) gives the 15 generators of an S0(6) subgroup, and 
likewise (Fy, ...., Fio) give the 6 generators of an S0(4) subgroup. These 
reahse the ^0(6) x 50(4), otherwise SU{4:), x SU{2)l x SU{2)r, Pati-Salam 
subgroup of SO(IO), the 3 + 3 generators of SU(2)l x SU(2)r being linear 
combinations of the 6 generators of SO (4). In terms of the creation and 
annihilation operators (j = 1,2,3) can give S0(6) and (j = 4,5) can give 
SO (4). We are concerned with transition operators X, a subset of those 
above, which form SO{<o) x SO{A) (equivalently SU{A)^ x SU{2)l x SU{2)r) 
invariant elements 

X 

by coupling to the scalar field (px which transforms in the same SU{A)f. x 
SU(2)l X SU(2)ji representation as {A\X\B). The neutral, colour singlet, 
subset Xq of X, coupling to vevs (px^ and multiplied by Yukawa coupling 
constants Yab{,Xq) 

Y.{A\X,\B)ct>x,YUX,) (5) 

yields contributions to 3 by 3 mass matrices. Here p, being 10 or 120 or 126, 
denotes the SO(IO) representation to which Xq belongs; thus the Yukawas 
are those appropriate to an unbroken SO(IO) symmetry. This symmetry is 
subsequently broken by the mass terms. 

The Table shows the subsets Xq contributing to the fermion masses via 
{A\Xq\B). The left hand column shows the dimensions of the S'f/(4)c x 
SU{2)i X SU{2)r representations to which the {A\ Xq \B) belongs. The right 
hand column gives the notation for the scalar vevs {(j)^) labelled also by the 
SO (10) representation (10, 120 or 126) of which the Xq is a member. For 
example the first pair of rows contain the neutral members of an SU{2)l x 
SU{2)ji bi-doublet and so do the second, third and fourth pairs of rows. 
Further comments on this Table as well as details of how the terms in it give 
rise to the mass matrices are given in the next Section. 
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Table 1: SO (10) Clifford algebra operators giving rise to fermion masses 



3 Yukawa coefficients, vevs, hermiticity 

The mass terms shown in the Table break not only the SO (10) symmetry 
but also the Pati-Salam left-right symmetry since the operators in the middle 
column do not commute with all the Clifford algebra Pati-Salam generators. 
As a theory of particle masses, in the present context, they are only appro- 
priate for use at high energy - say near the GUT energy - since the use of 
SO(IO) Yukawas as in equation ([5]) is only appropriate near energies where 
SO(IO) is a good symmetry. Such masses can ordinarily only be derived from 
experimental observation by the use of RGE equations such as those we shall 
make use of. 

The symmetry breaking is fairly clear. For example the last row can not 
only supply the heavy Majorana neutrino mass necessary for the Type I see- 
saw but also breaks the SU{2)ji symmetry because of the triplet component. 
It may be noted that the ninth row, necessary for the Type II seesaw breaks 
SU{2)i but because of the numerics of the neutrino masses only by a tiny 
amount. The combination (rir2 + r3r4 + T^Tq) occuring in rows 5,6,7,8 is 
proportional to the B-L generator and thus breaks the SU{A)c symmetry: 
SU{4)c -^ SU{3)c X U{1)b~l- In addition the operators in rows 9 and 10 do 
not commute with {B — L) so a non zero value of either of the corresponding 
vevs gives SU{4:)c — > SU{3)c- 

Thus if all the vevs in the Table are non- zero many symmetries are mul- 



tiply broken. 

The lines of the Table with one, three or five F^ correspond respectively 
to subsets of the {10}, {120} or {126} algebras of SO(IO). Fixing conventions 
we illustrate by outlining the down quarks mass matrix calculation. 

3.1 One r 

. From the {10} of SO(IO) these are the first two lines of the Table pro- 
viding the colour singlet members of an SU{2)l x SU{2)r bi-doublet. To 
corespond to colour singlet, neutral vevs Xq must be some combination of Fg 
and FiQ. The combinations Xs and Xs (being (Fio ±iF9)/2) are chosen with 
the corresponding vevs, v}^ = ((0io ± ^^g))- 

Evaluation of the resulting expressions {Y\^ being Yukawa coefficients) 

r+^\ BrC-\x5v'_' + xlvl') \^+b) Yl'j, + h.c. (6) 

gives directly contributions to the mass matrix elements of the quarks, charged 
leptons and neutrinos as can be seen on inspection of eqns. ([2]), ([3]) and (jl]), 
using the operator algebra. Consider for example the contribution to the 
two down quark mass matrix elements A,B and B,A. For one F these are all 
included in 



r-,B\ BvC-\, \^1,+a) YZ]v'^ + h.c. (7) 

This is expressed in terms of left, d^, and right, (ifi, Weyl spinors using C~^ = 
0"2 and equations ([2]), ([3]) with the correspondences ipa = cr2{daR)*]ipa5 = dai, 
a being the colour index. The result is 

{d^j^dsL + dij,dALy_\Yi'^ + yZ) (8) 

and the hermitian conjugate adds 

{di^dAR + d\^dBR) [v^'iYll + YZ)r (9) 

Obviously Yab = Yba follows without loss of generality. These two equations 
can be combined and writen in terms of 4-component Dirac spinors, dA and 
dR as 

2[Re{v^^YXl)dAdB + iIm{v^^YX%)dAl^dB\, (10) 



2[Re{v]_'Yll)dBdA + tIm{v]_'Yll)dBl5dA]. (H) 

and if v}PY}^ is real the pseudoscalar term vanishes and there is a real Dirac 
mass contribution to a symmetric flavour mass matrix. The other off-diagonal 
elements are of course completely similar and the diagonal contributions are 
real and scalar if f i°Y^|4 are real. Thus the flavour mass matrix is Hermitian. 
Also the mass contribution resulting conserves CP as well as P. 

3.2 Three T 

These give the {120} anti-symmetric representations which are developed 
here as imaginary hermitian matrix representations. 

(i) As shown in the Table there are two types of colour singlet neutral 
vevs associated with three T's. In what follows the evaluation for the type 

is outlined. The flavour mass matrix elements for the down quarks are 
included in 

[r^^l BrC-'TjTsX^ |^+^) Y^'^]vl''^ + h.c. (12) 

corresponding to equation ([7]). The total result of the calculation is 

- ^[(4i.rfBL - d^BRdAL)y_''''{Yr^ - Fi^°) 
HidLdAR - d^ALdsB)] [vi^'^Yll' - Y^'^)Y (13) 

Obviously Y^^ = —Y}^^ follows without loss of generality. These two 
lines can be rearranged and written in terms of 4-component Dirac spinors, 
dA and d^ as 

- 1^Re(y^:^''''Yf^)dAdB + iIm{v^}'^'^Yfi)dAl,dB\ (14) 

+ 2z[i?e(t;i2°"rif )JBrf^ + iIm[v^:^'^^Yf^)dsi^dA\ (15) 

Suppose v^]:^°'Yab fo be real. Then there is only a scalar mass term, 
the same applying to all the other off-diagonal contributions. The diagonal 
terms are anyway zero and thus the total contribution is imaginary anti- 
symmetric, so of hermitian flavour matrix form. 

The extra factor % with a SF operator such as FyFaXs is because FyFg = 
^(X4X4 ~ X4X4)- Thus it might be thought that the association of hermiticity 
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with a purely scalar mass term is due to a particular choice of phase in 
ryFgXs. This is not so; multiplication of that operator by e*"^ still yields the 
same association This can be shown by explicit calculation but generally one 
can reason as follows. 

The mass term in the Lagrangian of 3-fiavoured chiral theories is 

iP^M^l + Am^^r. (16) 

where M is a 3 x 3 flavour matrix and ^ij, tpi are 2-component spinors with 
the three flavour index implicit. For Hermitian matrices, M^ = M, this 
converts trivially into 4-component spinor, tp, as a scalar, ^M^, and upon 
flavour diagonalization results in normal Dirac equations. However if M is 
not Hermitian 

iPiM^L + AmHr. = ^^{M + M^)tlj - ^^75(M - Mt)V' (17) 

and the pseudoscalar parity breaking 75 term intrudes upon the canonical 
Dirac equation. Conversely the presence of a 75 term breaks hermiticity. 

(ii) Using the notation TV = Tir2 + r3r4 + FsFg it is seen from the Table 
that the other type of 3F operator is 

rFx5,rFxl 

For the down quarks this again yields the expressions flT^ - flT5]) with v^^'^ 
instead oiv]^^"'. Thus the association of hermiticity with a purely scalar mass 
term holds here too; and it equally applies for the masses of all particles. 

For three F in both(i) and (ii) above and contrary to the cases of one F 
and five F the adopted purely scalar case violates CP because of the anti- 
symmetry of the {120}. 

3.3 Five T 

These arise from the {126}. As shown in the Table mass terms are generated 
by 

TrrsFFxs and FyFgFFxs 

These indeed give results like the one F case in the respect that if v^^"'Y)]^ 
is real the pseudoscalar term vanishes and there is a real scalar mass con- 
tribution to a symmetric flavour mass matrix. That is a scalar mass matrix 
contribution implies that it is hermitian, and vice versa. As before the same 
applies to all the particle masses given by these operators. As in the One F 
case the resulting contribution also conserves CP. 



3.4 summary 

Collecting now the contributions from just the first eight lines of the Table 
these also include the familiar operator coupling right chiral neutrinos to left 
chiral neutrinos which is here denoted by M". All these together result in 
the mass operators with obvious labels 

Mis = 2Yllvl' - GYll'v'f^ + 2^ri|%{_20a + Q^Yll'v'f' (19) 

M^s = 2Y}lvl^ - 2rjf t;i26« _ 22FJ|%i^0" - 2tY}^^v'^^' (20) 

where the Yukawas Y^^ and Y^^ are real and symmetric while the Y\^ 
are real and anti-symmetric. The vacuum expectation values v are also real 
to make all the above mass matrices hermitian. 

It should be noted that the proof of hermiticity only evidently carries 
through in the fermion basis implicit in the present model; chiral rotations 
different for different flavours can change the coefficients of the scalar and 
pseudoscalar bilinears, changing the flavour mass matrices. 

The last two SF lines of the Table do not contribute mass to the quarks 
or charged leptons, but only to the neutrinos. This arises from evaluation 
of equations ([3]) using ([1]) and ([2]) and the associated particle assignments 
therein. These are given in the Appendix and associated with the left-right 
neutrinos of eqn 019p . First order block diagonalization of the resulting 6x6 
flavour mass matrix gives 3 light Majorana neutrinos and 3 heavy. 

4 Mass Relations 

The 16-plet of Weyl fermions eqns. (TTPj) has chiral neutrinos, n^, (n^jJ^.We 
denote the 4-component neutrinos corresponding to rii, ur as vl-, Nr re- 
spectively (see the Appendix for further details). The see-saw hypothesis 
assigns a large mass to Nji through a Majorana mass term, with flavour 
matrix (arising from eqn. ([3]) and the last line of the table) here denoted M. 
These left and right neutrinos couple together with flavour matrix M", given 
in eq. lflQl) . analogous to the quark and charged lepton mass matrices. The 
penultimate line of the Table gives rise to a Majorana mass term for vl The 
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flavour matrix for this, numerically very small compared to M, is denoted 
by m. Diagonalization of the resulting matrix 



Tijueutrinos 



m M" 
M"^ M 



(22) 



results in 3 light physical Majorana neutrinos and 3 heavy physical neutrinos 
in the top left and bottom right respectively of M"^"*''*"^'^. The light neutrino 
Lagrangian has: (i) a contribution having P and CP invariance containing a 
3x3 Hermitian flavour mass matrix; (ii) a contribution violating both P and 
CP containing a 3 x 3 imaginary non-Hermitian flavour mass matrix. This 
latter term arises from the original ('Type I') seesaw mechanism and is the 
only P violating term at high energy in the resulting mass Lagrangian of the 
model. 

For analyzing, as now follows, the mass relations of the quarks and 
charged leptons the notation is simplified, dropping the flavour indices from 
the Yukawa matrices, Y, and defining new matrices which incorporate the 
scalar vacuum expectation values (such as v}^ for those in the {10}). 

h = Y^%]_\ f = rl26a^m^ g ^ yl20^ma_ (23) 

Certain ratios of the real vacuum expectation values, required for the mass 
equations fl2Bl) to fl2^ below, are 

n = t;i206/^1^20a^ ^^ ^ v'f'/v'f\ (24) 

In addition for the Majorana neutrino matrices M = tmJ and m = Tmf 

ru = v'^^'/v^^^^ rm = v'^'^'/vl^'- (25) 

With these notations the quark and lepton mass matrix equations are 

M'^ = h-f -ig{l + ri) (26) 

M" = h + 3f - ig{l - 3ri) (27) 

M'' = rhh + rff + irgg{l-r2) (28) 

M"" = rhh-3rff + irgg{l + 3r2) (29) 

where the mass matrices are hermitian with h, f being real symmetric and g 
real antisymetric. As noted in the introduction the hermiticity is a significant 
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difference from the majority of previous papers. (Allowing for changes due to 
conventions the part of these equations in h and / are recognisably the same 
as those written in very many previous papers such as references [2] - [8] . On 
the other hand the terms in g are formally, and physically, different as they 
involve the ratios ri and r2 arising from the equations of section III.D.) 

Putting aside neutrino masses and mixing to be considered later, the 
present data, some of it being significantly only approximate, is 9 quark and 
charged lepton masses, 3 CKM matrix angles and 1 phase. 

To make use of the CKM matrix [V = UlUu) data it is a common de- 
vice, when possible, to take a basis in which either the d-quark matrix is real 
diagonal (implying Ud is unity) or the u-quark matrix is (so Uu is unity). 
Then either Uu or f/J respectively is the CKM matrix. While this might be 
done in the general case of the present model the unitary matrix required 
to change the basis generally bestows imaginary parts on the real matrices 
h,f,g, thus upsetting a simplifying feature of the model. However we can 
avoid this upset and shall make use of this device by considering some special 
cases. 

4.1 Special Cases 

To retain those features and the (relative) simplicity of numerical calculations 
there are two special cases of the vevs associated with the two couplings 
(which we have denoted as 120a and 120b respectively) in the {120}. These 
are (i) v™^ = y^^oa _^ ^^ _ -^^ making M" real symmetric, diagonalisable by 
a real orthogonal change of basis; and (ii) v}^^'^ = — ti^^oa ^ ^_^ _ _i^ making 
M'^ real symmetric, diagonalisable by a real orthogonal change of basis. The 
latter special case involves importing a relative phase of vr which however 
preserves the reality conditions of the model. Taking either of these special 
cases, with their associated change of basis, preserves the real symmetric 
(anti-symmetric) nature of the matrices h, /, (g). 

To illustrate the numerical evaluation consider the special case (ii) which 
will be used in the following section on neutrino masses and mixing. The 
mass matrix equations simplify to 

M'^ = h- f (30) 

M^ = h + 3f- Aig (31) 

M'^ = rhh + Tff + irgg{l - r^) (32) 

M" = Thh - Srjf + irgg{l + 3r2) (33) 
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The first three equations yield the mass relation 

M" = xM'^ + yReiM"") + izIm{M'') 



(34) 



which,equating matrix coefficients in the real and imaginary parts of equa- 
tions (13nD,(lSID, (IS2D, has 



X = {rj -3rh)/irh + rf), 



z 



-4/rg(l-r2) 



(35) 



The three matrices on the right hand side of equation flM]) are evaluated 
as follows in terms of the quark masses called here d, s, b and u, c, t,and the 
CKM matrix V. Now M'^,Re{M'^),Im{M'^) are written respectively as 



d 0' 




" Wi 


Ml 


Uz ' 







Vl 


vz 


s 


; 


Ml 


W2 


U2 


•) 


-Ml 





V2 


6 




Us 


U2 


W3 




-Ms 


-M2 






As stated above since the basis is diagonal in the down quarks and all the 
matrices are Hermitian then M" is diagonalized by the CKM matrix, V, 
giving the real parameters Ui,Vi,Wi in terms of u,c,t and the CKM matrix 
by 

M" = V^M:^,^^V (36) 

the diagonal matrix M^^^g having elements m, c, t. 

The choice of the CKM matrix is the additional input at this stage. 

Given a matrix, (i) its trace, (ii)the trace of the inverse and (iii) the deter- 
minant are each invariant under a unitary transformation. This leads to the 
known way of solving equations such as fl5^ since a unitary transformation 
can transform the Hermitian M*^ into a diagonal matrix of the charged lepton 
masses. Then we equate the three invariants of the two sides of equation (!M|) . 
For the right hand side (i) gives an expression linear in x, y while (ii) being 
equivalent to the sum of the 2x2 diagonal sub-matrices gives a quadratic 
in (x, y, z) but with z only occuring as z"^ and (iii) yields a cubic but again 
with z only occuring as z^ . 

The coefficients are real and elimination gives a cubic equation, either 
in X or y, which can be solved by an analytic expression and on evaluation 
yields numbers for the set x, y, z'^. Only those sets with x, y real and z^ > 
are acceptable solutions of equation ([31]). 
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5 Neutrino Masses and Mixings 

Continuing with the special case where ri = — 1 the hght neutrino masses 
and mixing are due for consideration. The approach is first to identify the 
three remaining free real parameters of the model. Then to choose a possible 
set of 3 neutrino masses which are in accord with mass data and search for a 
solution of those parameters (in principle there could be more than one) for 
which the model gives the chosen masses. This can be repeated for various 
choices of light neutrino masses. For each solution set the model is fixed with 
all parameters determined by the assumed quark, charged lepton and light 
neutrino masses and CKM matrix. Each parameter-fixed model predicts an 
MNS mixing matrix which can be considered for plausibility. 

So far use has been made of 13 assumed data points - 9 quark and charged 
lepton masses, 3 CKM angles and 1 phase - albeit some of these are subject 
to considerable error. Using equations (!30l) . (!3T|) . (!32|) . (!36!) and a basis diag- 
onal in the mass matrix M'^ the hermitian quark and charged lepton mass 
matrices M'^,M'^,M'^ have been synthesized from the 13 data points. An 
alternative description, in terms of the quantities on the RHS of equations 
f l30|) . fl3Tl) . fl32l) . fl33|) is that the real matrices h, /, g along with the real parame- 
ters r/i, Vf, rg{l — r2) have been constructed. (There seem to be 15 constructed 
parameters from 13 data points. However these 15 are not independent be- 
cause eq.( 15^ . arising from the special assumptions, produces 9 parameters 
from 7 data.) 

Turning now to the neutrinos, a httle manipulation of the equations yields 
the neutrino mass matrix of eq. (!33|) as 

iVf^ = ^—^M'' - ^-^4M + zl±^/m(M") (37) 

y y 1-^2 

M = {Re{M') - M'^)/4 = / = M/tm = m/r^. (38) 

In eq.( !37|l r2 occurs. It is a free parameter since only rg{l — r2) has been 
constructed. 

The 'seesaw' mass matrices m and M of eq. (l22l) (arising from the {126} 
representations of the last two rows of the Table) are proportional to the 
matrix / as displayed in fl38l) . Thus there are 2 more free parameters, r^ 
and tm, in addition to r2, leaving the model with just 3 so far undetermined 
real (and dimensionless) parameters. These can in priciple be fixed by fitting 
3 given light neutrino masses. Since data is only known on mass differences 
exploration requires postulation of one light neutrino mass. 

14 



By a redefinition of tlie neutrino states, working to first order in the 
numerically small matrix t] = M^M~^ ^ as shown in the Appendix, ]\^"eMtrj?ios 
can be transformed to block diagonal form 

rriy 
M 

The 3x3 matrix of the small mass neutrinos is given by the usual seesaw 
formula (Type II and Type I) rrii, = m—M^M~^M"''^. For the present special 
case the matrix m = r^^f, contributed wholly by the {126}, is hermitian, 
being real and symmetric. The Type I contribution, —M"'M~^M"''^ , though 
composed from hermitian matrices is neither real nor hermitian. However 
it is symmetric thus m,ycan be written in terms of real symmetric matrices 
■nil, ^^2 as rriu = mi + im2. The eigenvalues of the hermitian matrix 

m^ml = (rrii + im2){mi — im2) (39) 

are the neutrino masses squared. 

The questions now are: (i) are the 3 free parameters sufficient to fit some 
possible sets of experimental neutrino masses; (ii) if so is agreement also 
found with our knowledge of the MNS matrix 

5.1 An Example 

The Yukawa matrices, F ^°, y^^^, y^^*' have been assumed throughout to be 
appropriated to the 10,126,120 SO(IO) representations respectively. While 
this paper is not committed to a GUT model nevertheless the quark and 
lepton masses used should be those appropriate to a high energy where par- 
ity and SO(IO) are restored, thus involving extrapolation by renormalization 
group equations. There are well known (and well used) extrapolations by 
Das and Parida [T^ including renormalization group equations of the stan- 
dard (nonSUSY) model, the 2 Higgs doublet model and the minimum super- 
symmetric model. These have the general feature of such extrapolations of 
shifting the quark masses by considerably more than the relative shift of the 
lepton masses. For the light neutrino masses the present low energy data is 
here used as a guide. 

Possible sets of extrapolated quark and charged lepton masses and CKM 
matrix using the results of Das and Parida [H] have been constructed and 
used in various papers. We make use of some of those previous works by 

15 



adopting the data sets (appropriate to an energy of 2 x 10^^ GeV of (i) Goh 
et al. [2] and (ii) Bertolini et al.[7]). However it is not at all the purpose 
in this paper to make an assiduous search for solutions plausible on some 
criteria. Rather it is to sample sparsely to illustrate some possibilties albeit 
in a special case (ri = —1) of the original model. 

For the present experimental data the 3-neutrino mixing scheme reviewed 
by B. Kayser [13] is used. The neutrinos being named as 1,2,3 the following 
central values of the difference of squared masses are adopted 



(Am^)2i = 8.0 X 10-^eV"^ {Am^)32 = 2.5 x 10-^eV% 

along with the assumption of a neutrino hierarchy mli < m^2 < "^^3 so that 

^12 = ^li + (A"^')2i (40) 

ml3 = ml^ + (Am2)32. (41) 

The computing code calculates theoretical values of these 3 neutrino masses 
squared: /i^(r2, r^, tm ), each of these three being functions of the free param- 
eters r2, Tm, tm- It is necessary to find values of (r2, rm,rM) so that the fif are 
equal to the m^j, {i = 1, 2, 3). This is done in an obvious way by inventing 
and optimizing a function whose extreme value (say zero) is reached when 
IJ,^ = mli for all (i = 1, 2, 3), thus achieving an exact solution. The principle 
of this procedure is not the same as finding the 'best fit' (for example by 
X^ minimization) of many parameters. Because of the complication of the 
numerical calculation of the neutrino masses it is almost necessary to use a 
non- derivative method. So a non-derivative simplex method [T5], operating 
in the three-dimensional space of (r2, r^, ta/ ) is used to find any zeros of the 
chosen function. The computer search is specialised to those regions of small 
Tm and large tm suitable to produce small neutrino masses. 

(i) The paper of Goh et al.[2] uses the Das and Parida extrapolation of 
masses to 2 x lO^^ GeV by the MSSM RGB {tan{/3) = 10 in Table II of [H]) 
and also quotes the real part of the CKM m matrix. These are used as data 
input in the present example. Using the Wolfenstein parametrization gives 
Im(CKM) as 

-.00320 " 

-.00059 .00074 , (42) 

-.00328 

In all this mli is a quantity of choice being part of any postulate on 
neutrino masses. A first idea is that it should be small but not negligable, 
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say of the order of 1.0 x lO^^eV"^. In this region the optimization code 
gives the desired equahty to very great accuracy. For ml-^ = 2.0 x lO'^eV"^, 
m^2 — 1-0 ^ 10~^eV^^, m^3 = 2.6 x 10~^el^^ then then the dimensionless free 
parameters of equations fl371) and fl38|l have the values 



^2 



d.id, Tfji 



2.1 X 10" 



,n./ 



5.9 X 10^1 



(43) 



Having fixed the free parameters, the only extra input having been the 
partly hypothesized neutrino masses, it is then possible to calculate the MSN 
matrix. Rather surprisingly one finds good agreement with presently ac- 
cepted features of ths matrix. 

The charged lepton mass matrix, M^, (now calculable) is hermitian in 
the model and so diagonalizable by a unitary transformation 

XlM'^X, 



diag 



(44) 

= nil + ^^^2, and the matrix for the 
is hermitian and can be diagonalized 



The light neutrino mass matrix, m^ 
masses squared (mi + im2){rni — i'm2 
by a unitary transformation as 

Xl(mi + im2){mi-im2)X^. (45) 

The diagonalization of my requires a bi-unitary transformation as 

Xl{mi+im2)Yy. (46) 

The MNS matrix, similarly to the CKM matrix, is defined as 

Umns = XlXy. (47) 

and U can be calculated in the present model using the values of the 3 free 
parameters fixed by using values of the neutrino masses as discussed above. 
In this particular case the result for the matrix of moduli squared of the 
elements of U {\Uev\ ) is 



(48) 



bearing a distinct resemblance to the postulated 'ideal' structure of this ma- 
trix in tri-bimaximal mixing 



.638 


.344 


.017 


.260 


.331 


.409 


.102 


.325 


.573 



2/3 1/3 
1/6 1/3 1/2 
1/6 1/3 1/2 



(49) 
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It may be noted that the 13, 23 and 12 elements are clearly within the 
range of many analyses of the experiments ^U\. (It is interesting that in 
this solution the mass matrix elements contribution from Type II are much 
greater than those from Type I except for the (3,3) where the contributions 
are of the same order of magnitude.) Similar results hold for all values of 
mil between about 1.6 x 10~^ and 2.4 x 10~^ in eV'^. 

(ii) Bertolini et al. [7], making some use of the Das and Parida [H] ex- 
trapolation of masses to 2 x 10^^ GeV by the MSSM RGE {tan{/3) = 10) 
have given quark and charged lepton masses and 3 CKM angles plus the 
phase angle. Bertolini et al. [7] have revised the central values of the ex- 
trapolated masses of the lightest quarks to 0.55, 1.24, 21.7 MeV whereas Goh 
et al.[2], as used in (i), have 0.72,1.5,30.0 MeV. The extrapolated masses 
of the three heaviest quarks and the charged leptons remain the same as 
those in [2]. Some off-diagonal elements of the CKM matrix also display some 
non-trivial differences from those of Goh et al [2] used in (i). 

In the fitting, of the theory to the pseudo-physical neutrino masses, some 
samples of m^^ were taken in and around the same region as in (i) above. 
Any that succeeded in fitting yielded significantly different parameters from 
those in (i) resulting in spectacularly implausible MSN matrices. Further 
numerical investigation seemed to show that each one of the data changes 
mentioned above had influence on the results. That is results are sensitive to 
changes in the least well known data. It should be emphasized that neither 
in case (i) nor case (ii) was there attempted extensive investigation of very 
many neutrino spectra. 

6 Summary 

The model makes use of SO(IO) but with multiplets restricted in number 
and kind. 

Firstly the philosophy is to adopt the 16-plet of fermions as the particles 
we know and only consider couplings arising from 16 x 16 = 10 + 120 + 126. 
Secondly to restrict those couplings so that the resulting mass matrices are 
hermitian at high energy; an argument is given that within SO(IO) this can be 
formulated naturally; also the restoration at high energy of parity invariance 
of Lagrangian mass terms implies hermiticity. The scalar vacuum expection 
values are classified in the Pati-Salam subgroup of SO (10) using the Clifford 
algebra representation; this distinguishes two reahsations of the {120} asso- 



dated with different vacuum expection values. Assuming the values for the 
quark and charged lepton masses and the CKM angles and phase leaves just 
4 undetermined real dimensionless parameters in the theory. Fixing one of 
these parameters to a special value, thus simplifying the calculation and the 
range of results, leaves just 3 free parameters. To fix these the hierarchical 
hypothesis for light neutrino masses is adopted, together with specifying the 
lowest neutrino mass (at various values) to give the 3 light neutrino masses in 
accord with existing data on neutrino mass squared differences. The theory 
then predicts the MSN matrix. In the case of one set of quark and charged 
lepton masses and CKM matrix [2|,and a range of lowest neutrino masses, 
the MSN matrix is, rather surprisingly, in accord with the tri-bimaximal mix- 
ing suggested by the data [IB]. This result does does not hold for another 
data set [7J, for the same physical quantities, that was tried. Solutions thus 
appear sensitive to changes in the less well known physical quantities. 

Acknowledgments. I thank David Sutherland for many valuable dis- 
cussions and Colin Froggatt for comments on the manuscript. 
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APPENDIX 

In addition to the neutrino content of the standard model which has just 
left chiral (and massless) neutrinos, the 16-plet of SO(IO) has both left chiral, 
TT-i, and right chiral, ur, neutrinos. These occur in the 16-plet vector ([2]) as 
V's and —'00 respectively, being 2-coniponent Weyl spinors.. 

The matrix elements (jlj) give rise to three mass matrices for the neutrinos. 
Two of these are the self couplings of ni and tir respectively. Both of these 
are Majorana mass terms. The third couples ni and ur and gives rise to the 
mass term M"which appears in equations (!29|) and (133|) . The terms are as 
follows. 

(i) The self couplings of ur arise from X1X2X3X4X5 as in the last line of 
the table. They are 

{nARy''C-\nsRrv':''Yl'^ + hc 

where A, B are generation indices. These can be put into 4 component spinor 
notation by defining 4 component right chiral spinors 

NlR = {^,nAR) 
giving the (Majorana) mass terms as 



-N^j^NbrMab + he 

Mab = 2v'^''yX'^ (50) 

Since v^'^^ and Y^^ are real M is a real symmetric matrix. 

(ii) The self couplings of ni arise from X = X1X2X3X4X5 ^^ i^i the penul- 
timate line of the Table. Analogously to (i) above the resulting 2-component 
(Weyl) spinor results can be expressed in 4 component spinor terms by defin- 
ing 

^AL = i^AL, 0) 
giving the (Majorana) mass as 



-j^AL^BLniAB + he 



ruAB = 2v'_'''Yll' (51) 
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Since ti^^^^ and Y^^ are real m is a real symmetric matrix. 

(iii) The couplings of n/j to n^ arise from X as in the first 8 lines of the 
Table. In 4-component spinors these give mass terms 

TJIENbrM'Is + he 

where M" is the hermitian mass matrix of equation (1291) . 

These neutrino mass terms can be written in matrix form, using TTZM^N^ = 
'N^M'^^ul HH, as 



1 r 



VZ N- 



R 



m 



M 



N 



R 



+ hc. 



The matrix M is very large compared to M" and m. So, working to first 
order in M~^, the seesaw method can be implemented by defining 



L = ul + iM-^m)%N, R^ = Nr + {M-^^m)uL 



(52) 



noting that M and m are hermitian matrices. The neutrino mass terms 
become 



1 r 
2 



L R" 



m„ 







m,. 





M 

m + m^ 



R 



+ hc. 



-M'^M-^M''^ 



(53) 
(54) 



Eqn (l54l) gives the original seesaw term, often known as Type I seesaw, m 
is hermitian and ml is complex symmetric. Denoted {rrii,, M) are the flavour 
mass matrices of the (light,heavy) Majorana neutrinos with Majorana fields 



(55) 



u = L + L'',N = R'' + R. 
with the mass terms for the light neutrinos being 

-Vm,y{l + 75)2^ + he 
= Vmv + T'[Re{my) + i'y^Im{mi,)]i'. (56) 

In eqn fl56|) the two first terms, containing hermitian matrices, are P and 
CP invariant but the ^75 symmetric matrix term violates CP as well as P. 
Thus P and CP violation in the light Majorana neutrino mass terms arise 
solely from the original seesaw mechanism involving the non-hermitian prod- 
uct of 3 hermitian matrices. 
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